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We present numerical results on the transmission proper- 
ties of the left-handed materials (LHM) and split-ring res- 
onators (SRR). The simulation results are in qualitative 
agreement with experiments. The dependence of the trans- 
mission through LHM on the real and imaginary part of the 
electric permittivity of the metal, the length of the system, 
and the size of the unit cell are presented. We also study the 
dependence of the resonance frequency of the array of SRR 
on the ring thickness, inner diameter, radial and azimuthal 
gap, as well as on the electrical permittivity of the board and 
the embedding medium, where SRR resides. Qualitatively 
good agreement with previously published analytical results 
is obtained. 

PACS numbers: 73.20.Mf,41.20.Jb,42.70Qs 



I. INTRODUCTION 

Very recently, a new area of the research, called left-handed 
materials (LHM) has been experimentally deruaastrated by 
Smith et al.u'EI based on the work of Pendry et al.LlU. LHM are 
by definition composites, whose properties are not determined 
by the fundamental physical properties of their constituents 
but by the shape and distribution of specific patterns included 
in them. Thus, for certain patterns and distribution, the mea- 
sured effective permittivity tcff and the effective permeability 
/Xeff can be made to be less than zero. In such materials, the 
phase and group velocity of an electro-magnetic wave prop- 
agate in oBposite directions giving rise to a number of novel 
propertiegj. This behavior has-ibeen called "left-handedness" , 
a term first used by Veselagca over thirty years ago, to de- 
scribe the fact that the electric field, magnetic intensity and 
propagation vector are related by a left-handed rule. 

By combining a 2D array of split-ting resonators (SRRs) 
with a 2D array of wires. Smith et al.D demonstrated for the 
firsjt time the existence of left-handed materials. Pendry et 
al.u has suggested that an array of SRRs give an effective 
/ieff, which can be negative close to its resonance frequency. 
It is also well knownfl'u that an array of metallic wires be- 
haves like a high-pass filter, which means that the effective 
dielectric constant is negative at low frequencies. Recently, 
Shelby et al.u demonstrated experimentally that the index of 
refraction n is negative fot.a LHM. Negative refraction index 
was obtaiaed analyticallyo-and also from numerically simu- 
lated datatj. Also, Pendrylij has suggested that a LHM with 
negative n can make a perfect lens. 

Specific properties of LHM makes them interesting for 
physical and technological applications. While experimental 
preparation of the LHM structures is rather difficult, espe- 
cially when isotropic structures are required, numerical simu- 
lations could predict how the transmission properties depends 
on various structural parameters of the system. It will be 



extremely difficult, if not impossible, to predict the transmis- 
sion properties of such materials analytically. Mutual electro- 
magnetic interaction of neighboring SRRs and wires makes 
the problem even more difficult. Numerical simulations of var- 
ious configurations of SRRs and of LHMs could be therefore 
very useful in searching of the direction of the technological 
development. 

In this paper, we present systematic numerical results for 
the transmission properties of LHMs and SRRs. An improved 
version of the transfer-matrix method (TMM) is used. Trans- 
fer matrix was applied to problems of the transmission of the 
electro-magnetic (ravl)|-waves through non-homogeneous me- 
dia many years agot3~Ej. It was also used in numerical sim- 
ulations of the photonic band gap materials (for references 
see Ref. 15). TMM enables us to find a transmission and 
a reflection matrices from which the transmission, reflection 
and absorption could be obtained. The original numerical al- 
gorithm was described in Ref. 13. In our program we use a 
different algorithm which was originally developed for the aaJ- 
culation of the electronic conductance of disordered solidaij. 

The paper is organized as follows; In Section^ we describe 
briefly the structure. We concentrate on the structure dis- 
played in Figure 1. In Section |l^ we present and discuss our 
results. The dependence of the transmission of the LHM and 
SRR on the electrical permittivity o f the metallic co mpone nts 



of our structure is given in Section HI A. In Section IIIE we 
present the dependence of the transmission of the LHM on 
the size o f the unit cell and size of the metallic wires. In Sec- 
tion [II C we show the dependence of the resonanc e freq uency 

Section [II D deals 



of SRR on the parameters of the SRR. 
with the dependence of the resonance frequency on the per- 
mittivity of the board and embedding media. In Section ^ 
we summarize our results and give some conclusions. Finally, 
in the Appendix A we give detailed description of the transfer 
matrix method. 



II. STRUCTURE OF THE LHM 
META-MATERIAL 

Both in the experiment and in the numerical simulations, 
the left-handed meta-materials consist from an array of unit 
cells, each containing one SRR and one wire. Figure la shows 
a realization of the unit cell that we have simulated. The size 
of the unit cell Lx x Ly x Lz and the size of SRR itself are 
of order of mm. Waves propagate along the z-direction. The 
SRR lies in the yz plane, and the wire is parallel to the y axis. 

As we are interested mostly in the transmission properties 
of the left-handed meta-material, the configuration as pre- 
sented in Figure la, should be considered as one-dimensional. 
Indeed, such meta-materials possesses the left-handed prop- 
erties only for the electro-magnetic wave incoming in the z- 
direction and even then only for a given polarization. Two- 
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dimensional structures have been realized in experimentg3i3, 
in which two SRRs have been positioned in each unit cell in 
two perpendicular planes. For such structures, left-handed 
transmission properties have been observed for waves coming 
from any direction in the xz plane. No three - dimensional 
structure has been realized so far. 

Figure lb shows a single square SRR ©f the type used for 
our simulations and also for experimentscl. The structure of 
the SRR is defined by the following parameters: the ring 
thickness c, the radial gap d, the azimuthal gap g and the 
inner diameter r. The size of the SRR is 



27r7 = 



2L. 



: 4c + 2d + r. 



(1) 



Another parameter is the thickness of the SRR itself (in 
the i-direction). This thickness is very small in the experi- 
ments 0.02 mm). We can not simulate such thin struc- 
tures yet. In numerical simulations, we divide the unit cell 
into X Ny X mesh points. For homogeneous discretiza- 
tion, used throughout this paper, the discretization defines 
the minimum unit length 5 = L^/N^- All length parameters 
are then given as integer of 5. This holds also for the thick- 
ness of the SRR. Generally, the thickness of SRR used in our 
simulations is 0.25-0.33 mm. Although we do not expect that 
the thickness will considerably influence the electro-magnetic 
properties of the SRR, it still could cause small quantitative 
difference between our data and the experimental results. 



III. STRUCTURAL PARAMETERS 

A. Metallic permittivity 

The existence |-,|-, of 

LIfM has been experimentally demonstratedoQ for structures 
that have resonance frequencies in the GHz region. In this 
frequency region, we do not know the exact values of electri- 
cal permittivity em of the metal. We know that Im em is very 
large, a«ii/or the Re em is large but negative. In our previous 
studiesE-J we have found that the resonance frequency vq of 
the LHM depends only on the absolute value of em. In fact vq 
reaches the saturated value provided that |em| > 10*. Since 
we do not know the exact values of the metallic permittivity, 
we have studied the transmission of the LHM with different 
values of em. In the results presented in Figure 2[-yve choose 
em = 1 + j Im em with diflferent values of Im em.t3 The last 
is proportional to a{uj)/ux^. For simplicity, we neglect the 
oj-dependence of Im em and consider Im em = 8000, 18000 
and 38000 for the three cases presented in Figure 2. For each 
case, we present results of transmission for different number 
(1 to 10) of unit cells. Notice that the higher imaginary part 
of the metal the higher is the transmission. Also the losses 
due to the absorption are smaller, as can be seen from the 
decrease of the transmission peak as the length of the system 
increases. This Msult is consistent with the formula presented 
by Pendry et al.u for the effective permeability of the system 



and the resonance frequency 
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where p is the resistance of the metal, is the size of the 
system along the x axis, cught is a the velocity of light in vac- 
uum and parameters r, c and d characterize the structure of 
SRR. They are defined in Figure lb. Notice that the damping 
term 7 — > as o" — > 00. Since the Im em is proportional to ct, 
7 is inversely proportional to Im em. Our numerical results 
suggest that it is reasonable to expect that the LHM effect 
will be more pronounced in systems with higher conductivity. 

In Figure 3, we present the frequency dependence of the 
transmission for SRRs with the same parameters as those in 
Figure 2. Notice that the transmission is more pronounced 
as the length of the system is increased. Note also that the 
resonance gap becomes narrower when Im em increases. This 
is in agreement with Eq. . The frequency interval, in which 
the effective permeability is negative, becomes narrower when 
the damping factor 7 decreases. 

In Figure 4 we show the transmission through the LHM, in 
which the SRR are turned around their axis by 90 degrees. If 
we keep the same size of the unit cell as that of Figure 2, we 
do not obtain any LHM peak in the transmission, although 
there is a very well defined gap for the SRR alone. It seems 
that for this orientation of the SRR there is no overlap of the 
field of the wire with that of the SRR. The results shown in 
Figures 4 and 5 are therefore obtained with a reduced unit 
cell of 3.66 X 3.66 x 3.66 mm (the size of SRR is still 3x3 
mm). The LHM transmission peak is located close to the 
lower edge of the SRR gap (shown in Figure 5). This is in 
contrast to the results presented in Figures 2 and 3, where the 
LHM transmission peak is always located close to the upper 
edge of the SRR gap. Finally, the gap shown for the "turned" 
SRR shown in the Figure 5, is deeper and broader than the 
gap for the "up" SRR. 

In fact, for the SRR "up" structure, we found that the 
transmission in the gap is always of order of 10~^ — 10~*. We 
can explain this effect by non-zero transmission from the p to 
s polarized wave (and back). If the transmission t(j) — > s) and 
t(s — > p) > 0, then there is always the non-zero probability 
oc t(j> —> s)t{s p) for the p-polarized wave to switch into 
the s state, at the beginning of the sample, move throughout 
the sample as the s wave (for which neither wires nor SRR are 
interesting) , and in the last unit cell to switch back into the p- 
polarized state. This process contributes to the transmission 
probability T{p — > p) of the whole sample and determines 
the bottom level of the transmission gap for SRR. We indeed 
found that t{s p) ~ 10~* for the "up" SRR. In the "turned" 
SRR case, both t{p — > s) and p) should be zero due to 

the symmetry of the unit cellE3 Our data give t{p ^ s) ~ 
t{s ^ p) ~ 10"'' for the "turned" SRR which determines the 
decrease of the transmission in the gap below 10~^^. 



Mcff 



i/q + i-yv 



(2) 



with the damping factor 
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B. Dependence on the size of the unit cell and the 
width of the metallic wire. 



[II A 



the transmission peak 



As we discussed in Section 
for the LHM with cuts in the SRR in the horizontal direction 
appears only when the size of the unit cell is really small. The 
effect of the size of the unit cell was demonstrated already in 
Figures 2 and 3, where we compared the transmission for the 
"up" SRR and LHM of different size of the unit cell. Both the 
transmission gap for an array of SRRs and the transmission 
peak for LHM are broader for smaller unit cell. 

In Figure ^, we show the transmission for the LHM struc- 
ture with q unit cell of size 5 x 3.66 x 5 mm for the "turned" 
SRR. Evidently, there is no transmission peak for all the sys- 
tem lengths studied. The size of the unit cell must be reduced 
considerably to obtain a transmission peak. 

Figure ^ presents the transmission peak for various sizes of 
the unit cell. Resonance frequency decreases as the distance 
between the SRRs in the x direction decreases. This agree 
qualitatively (although not quantitatively) with theoretical 
formula given by Eqn. (^). We see also that an increase of 
the distance between SRR in the z direction while keeping 
the Lx constant causes sharp decrease and narrowness of the 
transmission peak. 



C. Resonance frequency of SRR 

In this section we study how the structure of the SRR influ- 
ences the position of the resonance gap. In order to simulate 
various forms of the SRR, we need to have as many as possible 
mesh points in the yz plane. Keeping in mind the increase 
of the computer time when the number of mesh points in- 
creases, we used a unit cell with Lx < Ly,Lz. The actual 
size of the unit cell in this section is ~ 2.63 mm and 
Ly — Lz — 6.05 mm. and we use uniform discretization with 
Nx X Ny X Nz = 10 X 23 X 23 mesh points. This discretization 
defines a minimum unit length S = 0.263 mm. SRR with size 
of « 5 X 5 mm, is divided into 19 x 19 mesh points. 

The electrical permittivity of the metallic components is 
chosen to be em — —1000 + 10.000 i. We expect that larger 
value of Im, em will increase a little the position of the reso- 
nance gaptj. However, the dependence on the different struc- 
tural parameters will remain the same. Higher values of |em|, 
however, will require more CPU time because of shorter inter- 
val between the normalization of the transmitted waves (see 
Appendix A for details). 

We have considered 23 different SRR structures and stud- 
ied how the resonance frequency fo depends on the structure 
parameters. The lowest vo — 3.75 GHz was found for SRR 
with c:d:r:g=l;l;13:l. On the other hand, SRR 
with c:d:r:g = 2;3:5:3 exhibits vq = 6.86 GHz. 

We present our results on the dependence of vo on the az- 
imuthal gap c (Figure radial gap g (Figure ^ and ring 
thickness d (Figure p^ ). In all these cases vg increases as the 
different parameters increase. The dependence shown in Fig- 
ures B-[ld agrees cpalitatively with those done by a different 
numerical methodtiJ. When compare our resultei with the an- 
alytical arguments presented by Pendry et alD, we have to 



keep in mind that that various assumptions about the struc- 
tural parameters have been done in derivation of Eqn. (^, 
which are not fulfilled for our structure. Note also that the 
azimuthal gap does not enter the formula for the resonance 
frequency given by Eq. (^. Moreover, as the size of the 
SRR is constant, the structural parameters are not indepen- 
dent each form other. Thus, due to the Eq. (|l|), increase of 
the azimuthal gap causes decrease of the inner diameter and 
vice versa. When taking these restrictions into account, the 
agreement with analytical results is satisfactory. 



D. Material parameters 

In Figures |ll]and0 we show how the resonance frequency 
depends on the permittivity of the dielectric board and on 
the permittivity of the embedding media. As expected, the 
resonance frequency decreases considerably with the increase 
of the value of both permittivities. 



IV. CONCLUSION 

In summary, we have used the transfer matrix method to 
calculate the transmission properties of the left-handed ma- 
terials and arrays of split ring resonators. The role of absorp- 
tion of the metallic components of our SRR and LHM has 
been simulated. It is found that the LHM transmission peak 
depends on the imaginary part of the metallic permittivity 
em, the length of the system and the size of the unit cell. 
Higher conductivity of the metal guarantees better transmis- 
sion properties of LHM. 

For an array of SRR, the resonance frequency i^o was com- 
puted and is found to agree with experimental data. The 
dependence of the resonance frequency on various struc- 
tural parameters of the SRR were numerically obtained and 
compared with analytical estimates and also with other nu- 
merical techniques. 

The main disadvantage of the presented transfer-matrix 
method is that it can not treat structures with smaller length 
scales than our discretization mesh. For example, the thick- 
ness of the SRR is an order of magnitude smaller in experi- 
ments than in our simulation. Also structural parameters of 
SRR can be changed only discontinuously as multiplies of the 
unit mesh length. This could be partially overcomed by gener- 
alizing the present code to a non-uniform mesh discretization. 
Nevertheless, already uniform discretization enables us to ob- 
tain credible data. Comparison of our resuJts with those ob- 
tained by the commercial software MAFIAE3 confirmed that 
both methods find the same position of the resonant gap pro- 
vided that they use the same mesh discretization. 

Our nuiBprical data agree qualitatively with the experimen- 
tal results.tJ As we can not tune the exact parameters of SRR 
(as well as its circular shape), and when taken into account the 
strong dependence of the resonance frequency on the permit- 
tivity of the board, we do not expect to obtain very accurate 
quantitative agreement with experimental data. 

Our studies demonstrate that the transfer matrix method 
can be reliable used to calculate the transmission and reflec- 
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tion properties of left-handed materials and split-ring res- 
onators. Thus, numerical simulations could answer some 
practical questions about different proposed structures, which 
might be too complicated to be treated by analytical studies. 
The transfer matrix method can be used in the future for de- 
tailed studies of two-dimensional and even three-dimensional 
structures. These structures should contain more SRRs and 
wires per unit cell, which makes their analytical analysis ex- 
tremely difficult. On the other hand, it is extremely important 
to find the best design and test the transmission properties 
of proposed meta-material even before their fabrication and 
experimental measurements start. 



APPENDIX A: 

Transfer matrix calculations are based on the scattering 
formalism. The sample is considered as the scatterer of an 
incoming wave. The wave normalized to the unit current is 
coming from the — oo, and is scattered by the sample. Scat- 
terer is characterized by four parameters: transmission of the 
wave from the left to the right {t+), from the right to the 
left (t-), and by reflection coefficient from the right to the 
right (r+) and from the left to the left (r_). Corresponding 
scattering matrix S reads: 



S = 



t+ r+ 
r- t- 



(Al) 



which determines the amplitudes of the outgoing waves B, C 
in terms of the amplitudes of the incoming waves A,D: 



Q 



Relation ( A2 ) can be re- written into the form 



(A2) 



(A3) 



T is the transfer matrix, which determines the fields on one 
side of the sample with the fields on the another side. Its 
explicit form reads 



tZ -tZ r- 
r+tZ^ t+-r+tZ^r- 



(A4) 



Transfer matrix (TM) fulfills the composition law. If the 
sample consists from two subsystems, then the transfer ma- 
trix 7i2 of the whole sample can be calculated form transfer 
matrices of its subsystems as 



712 ~ T2T1. 



(A5) 



Resulting TM has again the form (A4). This composi- 
tion law enables us to calculate transmission of complicated 
structure from the transfer matrices of its parts (thin slices). 

In numerical calculations, the total volume of the system 
is divided into small cells and fields in each cell are coupled 
to those in the neighboring cell. We discretize the Maxwell 
equations following the method described in Ref. 13. In each 
point of the lattice we have to calculate four components of 

Hy. 



We assume that our system is connected to two semi- 
infinite leads (with e = 1 and jj, = 1). EM wave is coming from 
the right and is scattered by the sample. Resulting waves ei- 
ther continues to the left on the left side lead, or are traveling 
back to the right on the right side lead. Periodic boundary 
conditions in the directions perpendicular to the direction of 
the wave propagation are used. 

We decompose the system into n thin slices and define a 
TM for each of them. Explicit form of the TM for a thin slice 
is in Refs. 12 and 13. The EM field in the {k + l)th slice can 
be obtain from the fcth slice as 



c&fe+i = 7I-$fc 



(A6) 



with Tk being the transfer matrix corresponding to the fcth 
slice. The transfer matrix T of the whole sample reads 



T — TnTn-l ■ ■ ■ T2T1 



(A7) 



If there is A'^ mesh points in the slice, then the length of 
the vector $ is 4A'' (it contains 4 components of EM field in 
each point). 

Note that we are able to find the explicit form of the TM 
only in the real space representation. To obtain the transmis- 
sion, we have to transform the TM into the "wave" represen- 
tation, which is defined by the eigenvectors of the TM in the 
leads. Therefore, in the first step we have to diagonalize the 
TM in the leads. 

Each eigenvalue of the TM is two time degenerate because 
there are two polarizations p and s of the EM wave. More- 
over, if A is an eigenvalue, then is also an eigenvalue 
corresponding to the wave traveling in the opposite direc- 
tion. In general, the TM has some eigenvalues with modulus 
equal to 1: A = expifc. The corresponding eigenvectors rep- 
resent propagating waves. Others eigenvalues are of the form 
A = exp ±K. They correspond to the evanescent modes. For 
the frequency range which is interesting for the LHM studies, 
the TM has only one propagating mode. 

As the TM is not Hermitian matrix, we have to calculate 
left and right eigenvectors separately. From the eigenvectors 
we construct three matrices: The 2N x 4A'^ matrix _Ri contains 
in its columns 2A'^ right eigenvectors which correspond to the 
wave traveling to the left. Matrices L\ and L2 are 4A'^ x 2A^ 
matrices which contains in their rows the left eigenvectors 
corresponding to waves traveling to the left and to the right, 
respectively. 



The general expression of the TM given by Eo— ( A4 ) en- 
ables us to find the transmission matrix explicitlytj 



tZ = LiTRi 
and the refiection matrix from the relation 
r+tZ^ = L2TR1. 



(A8) 



(A9) 



At this point we have to distinguish between the propagating 
and the evanescent modes. For a frequency range of interest, 
the TM in leads has only one propagating mode for each di- 
rection. We need therefore only 2x2 sub- matrices t_ (ij) and 
r+{ij) with i, j = 1 or 2 for the p or s polarized wave. The 
transmission and reflection are then 



7,, 



t.{ij)t*_{ij) 



■ r+{ij)r*+{ij) 



(AlO) 
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and absorption 

Ap - 
As -- 



- T 



Rs 



Rst 



(All) 



It seems that relations (AJ) and (A9) solve our problem 



completely. However, the above algorithm must be modified. 
The reason is that the elements of the matrix tZ^ are given by 
their larger eigenvalues. We are, however, interesting in the 
largest eigenvalues of the matrix t_. As the elements of the 
transfer matrix incr ease exponentially in the iteration proce- 
dure given by Eq. (A7), an information about the smallest 
eigenvalues oi tZ^ will be quickly lost. We have therefore to 
introduce some re-normalization procedure. We use the pro- 
cedure described in Ref. 16. 



Relation (A.8) can be written as 



t~_ 



Lir 



(n) 



(A12) 



where we have defined 2N x 4A'' matrices r'*', k — 0,1, . . .n 
as 



^Ri 



(A13) 



Each matrix r can be written as 



^^^^^ 

with ri, r2 being the 2A'' x 2A'' matrices. We transform r as 

(A15) 



: r n 



1 

r2ri 



and define r'-'"-' = Tk{r')^''~^\ In contrast to ri and r2, all 
eigenvalues of the matrix r2rZ^ are of order of unity. Relation 
( A12 ) can be now re- written into the form 



Li 



J") 



„(")_("-!) 
' 1 ' 1 



from which we get that 



JO) 



1 



1 



(A16) 



(A17) 



From Eqn. ( A.£ ) we find 



/ O '1 



„(") 



(A18) 



The matrix inversion in the formulae ( A15 - A18 ) can ob- 
tained also by the soluton of a system linear equations. In- 
deed, matrix BA~^ equals to matrix X, which solves the sys- 
tem of linear equations B — XA. CPU time could be reduces 
considerably in this way, especially for large matrices. 



All elements of the matrices on the rhs of Eqn. (A17) are 
of order of unity. The price we have to pay for this stability is 
an increase of the CPU time. Fortunately, if the elements of 
the transfer matrix are not too large (which is not the case in 



systems studied in this paper), then it is enough to perform 
described normalization procedure only after every 6-8 steps. 
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D. Vier, private communication 




FIG. 1. Top: The structure of the unit cell as was used 
in the present simulations. Structure acts as the left-handed 
meta-material if the electro-magnetic wave propagates along 
the z direction and is polarized with electric field E parallel 
to the wire and magnetic field H parallel to the axis of SRR. 
Bottom: The structure of the SRR and definition of the SRR 
parameters. 
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FIG. 2. LHM peak for various values of the metallic per- 
mittivity. Lines corresponds to system length of 1,2,. . . 10 unit 
cells. SRR is modeled as in Figure |l| The size of the unit 
cell is 5 X 3.66 x 5 mm, the size of SRR is 3 mm, and the size 
of the wire is 1 x 1 mm. The dashed line is transmission for 
LHIVI system with unit cell 3.66 x 3.66 x 3.66 and the system 
length of 10 unit cells. 
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FIG. 3. Transmission for SRR for the same systems as in 
Figure ^ Data confirm tliat tlie resonance frequency d(f«p 
not depend on the metallic permittivity. This agrees withE-J. 
However, resonance gap becomes narrower as Im em of the 
metallic components increases. Shaded area represents a gap 
for the array of the SRR with a unit cell 3.66 x 3.66 x 3.66 
and length system of 10 unit cells. Note also that the depth 
of the transmission is constant (~ 10~^) and see text for ex- 
planation. 
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FIG. 4. Transmission for LHM with the SRR rotated by 
90 degrees. The size of the unit cell is now 3.66 x 3.66 x 3.66 
mm. We found no transmission peak for the unit cell size as 
in Figure ^ (see Figure 11). For comparison with LHM with 
"up" oriented SRR, we show also the peak for this structure 
of the same unit cell and metallic permittivity e = 1 + 38.000 i 
and length system of 10 unit cells (dashed line). 
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FIG. 5. Transmission for SRR rotated in 90 degrees. The 
size of the unit cell is now 3.66 x 3.66 x 3.66 mm. For com- 
parison, we add also data for the array of "up" SRRs and 
tin = 1 + 38.000 i and length system of 10 unit cells (shaded 
area) . 
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FIG. 6. The transmission for LHM for the same parameters 
as in Figure 4, but with larger unit cell 5 x 3.66 x 5 mm. No 
transmission peak is observed. 
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FIG. 7. Transmission peak for various sizes of the unit cell 
with "turned" SRR. The system length is 10 unit cells. Res- 
onance frequency decreases slowly as the distance of SRR in 
the a;-direction decreases. Increase of the distance of SRR 
along the z direction causes decrease of the transmission peak, 
which becomes also much narrower. 
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FIG. 8. Resonance frequency as a function of the azimuthal 
gap g for various SRR structure. The size of the SRR is 5 mm, 
the unit cell is 2.63 x 6.05 x 6.05 mm, which corresponds to 
10 X 23 X 23 mesh points. The ratio d : c : r is 1 : 1 : 13, 
3 : 1 : 5, 2 : 1 : 9, 1 : 2 : 11 and 2 : 2 : 7, respectively. Note 
that azimuthal gap does not enter in Eqn. (M). 
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FIG. 9. Resonance frequency as a function of the radial 
gap c (distance between rings) for various SRR structures. 
The size of the sample and of the unit cell is an in Fig- 
ure ^. In contrast to Eqn. (^), which predicts decrease of 
i/Q when c increases, we found increase of the resonance fre- 
quency. This could be explained by the fact that an increase 
of the radial gap causes decrease of the inner diameter because 
r + 2c = 3.947 mm (o) and 2.895 mm (A). Presented data 
can be therefore interpreted also as the inner diameter - de- 
pendence of the resonance frequency. The last {uq ~ r~^^^) is 
much stronger than the logarithmic dependence vq ~ ln~^'^^ c. 
Then, presented data confirm that the decrease of the inner 
diameter causes an increase of the resonance frequency. 
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FIG. 11. Resonance frequency as a function of the permit- 
tivity of board. The size of the sample and of the unit cell is 
an in Figure B. 
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FIG. 12. Resonance frequency as a function of the permit- 
tivity of embedding media in the unit cell (the permittivity 
of semi-infinite leads remains 1). The size of the sample and 
of the unit cell is an in Figure 
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FIG. 10. Resonance frequency as a function of ring thiclc- 
ness for two sets of SRR. The size of the sample and of the 
unit cell is an in Figure ^. 
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